ASYMPTOTIC ANALYSIS OF ORTHOGONAL POLYNOMIALS VIA THE 

TRANSFER MATRIX APPROACH 



MANWAH LILIAN WONG* 



Abstract. In this paper, we present a new method via the transfer matrix approach to obtain 
asymptotic formulae of orthogonal polynomials with asymptotically identical coefficients of bounded 
variation. We make use of the hyperbolicity of the recurrence matrices and employ Kooman's 
Theorem to diagonalize them simultaneously. The method introduced in this paper allows one to 
consider products of matrices such that entries of consecutive matrices are of bounded variation. 

Finally, we apply the asymptotic formulae obtained to solve the point mass problem on the real 
line when the measure is essentially supported on an interval. We prove that if a point mass is added 
to such a measure outside its essential support, then the perturbed recurrence coefficients will also 
be asymptotically identical with the same limit and of bounded variation. 

1. Introduction 
Let [i be a non-trivial measure on M. such that for all n G N, 

\x\ n dfi(x) < oo. (1.1) 
We form an inner product and a norm on L 2 (M, d[i) as follows: for any f,g £ L 2 (JH,dfi), define 

(/,<?> = / 7(x)9(x)d f i(x), (1.2) 
f(x) 2 dfi(x). (1.3) 



By the Gram-Schmidt process, we can orthogonalize l,x,x 2 ,... and obtain the family of or- 
thogonal polynomials. We denote the n-th monic orthogonal polynomial as P n (x) and the n-th 
orthonormal polynomial as p n (x). Let 

«n = TE~\\- (1-4) 
Then the n-th orthonormal polynomial is given by 

P n (x) 

Pn(x) = = n n x n + lower order terms. 0-5) 

1 1 *Tl 1 1 

It is well-known that these orthogonal polynomials satisfy the following three-term recurrence 
relations: 

xP n (x) = P n+ \(x) + b n+ i P n (x) + a 2 P n _i(x) (1.6) 
xp n (x) = a n+ i(x)p n+ i{x) + b n+1 p n (x) + a n p n -i(x) (1.7) 



Date: March 31, 2011. 

2010 Mathematics Subject Classification. 33D45, 42C05, 34E10. 

Key words and phrases, point mass, asymptotic analysis, orthogonal polynomials, bounded variation, Nevai class. 
* School of Mathematics, Georgia Institute of Technology, 686 Cherry Street, Atlanta, GA 30332, USA. 
(wongOmath . gatech . edu) . 

1 



2 



M.-W. L. WONG 



with the properties that 



\Pn\ 



\Pn- 



>0 and b n+ i = (xp n ,p n ). 



(1.8) 



The (a n , b n )^ =1 are called the recurrence coefficients associated to the measure d\x. The recurrence 
relation (|1.7p is often represented by the matrix 
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(1.9) 



which is called the Jacobi matrix. 

Apart from the Jacobi matrix, there is another representation of the recurrence relation by means 
of the transfer matrix T n (x): observe that the recurrence relation (jl.7p can be rewritten as 

a n+ ip n+ i(x) = (x- b n+1 )p n (x) - a n p n _i(x) for n > 0. (1.10) 

Therefore, the general solution of (jl.lOp above can be expressed in the following way: 

K a n p n (x) 



An+l(x) 



Pn\X) 



1.11 



where 



Aj(x) 



-bj -1 
a) 



(1.12) 



This motivates the definition of the transfer matrix 

T n {x) = A n (x)A n _ 1 (x) ■ ■ ■ Ai(x) for n > 1. 

In particular, the recurrence relation could be expressed in terms of the transfer matrix applied to 
(l,0f: 

Pn(x) 

t _ip n _i(x) 



T n {x) 



;i.i3) 



1.14) 



Moreover, since det Aj{x) = 1 for all j > 1, 



1.15) 



detT n (x) = JJdetAj(a:) 

The transfer matrix T n (x) and the ^4 n (x)'s will serve as important tools when we derive the 
asymptotic formulae for the p n (x) 's in Section [H 

The reader should be reminded that besides the definitions given for T n (x) and A n {x) in (11131) 
and (1112p . there are a few other commonly used definitions in the literature. 

Let J n be the truncated n x n matrix obtained from the first n rows and columns of the Jacobi 
matrix J. Note that by the recurrence relation, 

/ Po( x ) \ 



(J n ~ x) 



pi(x) 

Pn-2{x) 
\p n -l{x)) 



( 










-a n p n (x)J 



1.16) 



Every zero of p n (%) is an eigenvalue of J n and the orthogonal polynomials form the eigenfunction 
for the truncated Jacobi matrix. In other words, p n (x) is the characteristic polynomial of J n . More- 
over, it is known that the spectral theory of one-dimensional operators (eg., the one-dimensional 
discrete Schrodinger operator) has a lot in common with the classical theory of orthogonal poly- 
nomials on the real line. Many results were proven for the Jacobi operator and the Schrodinger 
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operator in parallel (see for example [5]). Therefore, orthogonal polynomials have gained the at- 
tention of both the spectral theory community and the Schrodinger operator community in recent 
years. 

For a more comprehensive introduction to the theory of orthogonal polynomials on the real line, 
the reader may refer to [H [HJ [T71 [23] . 

2. Results 

First, we consider a measure in the Nevai class At (a, 6), which consists of measures with re- 
currence coefficients satisfying a n —> a,b n — > b. It is well-known that measures in A4(a,b) have 
essential support [b — 2a, 6 + 2a] (see [T7J for a detailed discussion). In this paper, we limit ourselves 
to a n — > a (Dombrowski [6j showed that a Jacobi matrix with liminf \a n \ = has empty a.c. 
spectrum). For the importance of the Nevai class and references to the many investigations thereof, 
the reader may refer to [151 HHJ H7J HB] . 

Apart from the asymptotic formulae listed in Theorem 12.11 below, it is interesting to note the 
method developed in this paper by means of applying Kooman's Theorem to the transfer matrix (see 
Section H]). We make use of the fact that if xq is outside the essential support of the measure, then 
A u {xq) is hyperbolic for all large n. This fact allows us to apply Kooman's theorem to simultaneously 
diagonalize A n and deduce asymptotic formulae for p n (x) outside the essential support. Due to the 
length of the proof, a sketch of the proof is provided in Section 14.11 

The first result can be summarized as follows: 

Theorem 2.1. Let fj, be a measure in A4(a, b) with recurrence coefficients of bounded variation and 
xq G M\[b — 2a, b + 2a]. Then the asymptotic formulae for p n (xo) ar ^ o-s follows: 

(1) ju(xo) > if and only if given any e > 0, there exists a constant C e such that 

\Pn(x )\<C e (\- +e) n . (2.1) 

where |A _ | < 1 is the eigenvalue of A qo (xq) (see (|4.ip for the definition of A^). 

(2) (j,{xq) = if and only if for every n G N, p n ( x o) *s * n the form 



Pn(x )= M A I \ k m (2-2) 




where (k n ) n is a convergent sequence of bounded variation that varies according to initial 
conditions, and is the eigenvalue of the recurrence matrix Aj(x) whose norm is great 
than 1 (see (j!.12[) for the definition of the recurrence matrix Aj(x)) . The sequence will be 
computed explicitly in the proof in Section^ 

Then we apply the asymptotic formulae obtained to solve the point mass problem. We add a 
pure point xo G M. to /i to form the measure fi as follows: 

fl = fi + ^S Xo , 7 > 0. (2.3) 

In Theorem 12.21 below, we give formulae relating the orthogonal polynomials and the recurrence 
coefficients of fi and fi. Even though these formulae are known (see pZ]), the proofs are given below 
for the convenience of the reader. 

A note on notation. We shall denote objects associated to the measure jl with a ~ on top. 
For example, the n-th monic orthogonal polynomial with respect to the measure jl is denoted as 

Pn{x). 

Theorem 2.2. Let [i and jl be measures defined as in (12. 3D . Then the n-th monic orthogonal 
polynomials of fi and jl are related by the following formula: 

jP n (x )K n (x,x ) 



Pn{x) ; K ' 



Pn(x) 



1 + jK n (x ,x ) 



(2.4) 
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with 



K n \ 2 1 +jK n (xo,x ) 



and the reproducing kernel K n (x,y) is defined as 



(2.5) 



K n {x,y) = J2Pj( x )Pj(y)- ( 2 - 6 ) 

3=0 

Furthermore, the recurrence coefficients of n and ft are related as follows: 
(1) 

a n = a n J > (2.7) 

V tra 

where 

1 + ~fK n _ 1 (x ,x ) 
1 + jK n (xo,xo) 



(2) 



On+1 = O n+ l — ■; r 1 — — — . (2.9) 

1 + 7A„_i(x ,x ) l+jK n (xo,x ) 



In Section [7] we combine those asymptotic formulae with Theorem 12.21 to prove the following 
result: 

Theorem 2.3. Let \x be a non-trivial measure on R such that its recurrence coefficients satisfy 

a n ^ a ^0, b n ^b; (2.10) 

oo 

\a n+ i - a n \ + \b n+ i - b n \ < oo. (2.11) 

n=0 

The essential support of \i is [b — 2a, b + 2a] . If we add finitely many distinct pure points 
xx, . . . , Xk G R\[6 — 2a, b + 2a] to /i as follows 

k 

= M + X] 7 -?^ ' 7 -J > °' ( 2 ' 12 ) 
£/ien t/ie recurrence coefficients of fik satisfy ()2.10|) and (|2.1ip . 

3. ASYMPTOTIC ANALYSIS AND THE POINT MASS PROBLEM 

The point mass problem has a very long history (see the Introduction of [26] for details) and it 
has its physical significance. As noted in the Introduction, results were often proven for the Jacobi 
operator and the Schrodinger operator in parallel. As a result, the point mass problem has been 
investigated by both the orthogonal polynomials and the mathematical physics communities. 

The earliest work related to the point mass problem could be due to Wigner-von Neumann |24j . 
where they constructed a potential with an embedded eigenvalue. In 1946, Borg [1] proved a well- 
known result concerning the Sturm-Liouville problem, which implies that if two Sturm-Liouville 
operators have spectra differing by a finite number of eigenvalues, then their corresponding potential 
functions might not be the same. Later, GePfand-Levitan [7] showed that in order to recover the 
potential one also needs the norming constants, which correspond to the weights of pure points in 
the context of orthogonal polynomials. 

The point mass problem has been considered under various settings. In |25| Szwarc considered a 
measure with bounded support S C [0, +oo) of which the recurrence coefficients satisfy a n a~, 1 — > 1 
and a n /b n — > y/A for some A > 0. It was proven that if a point mass is added to the measure, then 
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a n — a' n — > and b n — b' n — > as n — > oo, where a' n , b' n are the perturbed coefficients, though the 
specific rate of convergence was not shown. 

Szwarc's paper also discusses the growth of orthogonal polynomials for measures supported on 
[0, 1] with a finite number of pure points, with a n —> 1 and b n — > 1/2. It was proven under such 
conditions, 

limsup|p n (x)| 1/n < 1 (3.1) 



n 



uniformly on [0,1]. Unfortunately, such results (mostly derived from |19|). like many existing 
asymptotic results on p n (x), are insufficient to prove Theorem 12.21 

In [16], Mate-Nevai-Totik considered a measure in the Nevai class with coefficients satisfying 
a n — > 1 , b n — > and 

\a n +i - a n \ + \b n+ i - b n \ < oo. (3.2) 

n 

This measure is supported on [—1, 1]. This is a special case among those considered in Theorem 12.31 
The authors gave asymptotic formulae for p n (x) on a compact set K C C with K D [— 1, 1] = as 
well as on a compact set K' C supp(/x)\[— 1, 1]. The authors also presented ratio asymptotic results 
concerning p n {x) / p n -\{x) . However, in both cases, we need more precise bounds on error terms to 
establish Theorem 12.31 (see Section 

There is an analog of the point mass problem on the unit circle, but the problem is different in 
nature as the recurrence coefficients of a non-trivial probability measure on <9B form a one-parameter 
family (a„)^L with a„eB (also known as the Verblunsky coefficients in [201 [21]). 

The problem on the unit circle reads as follows: let v be a non-trivial probability measure on 
<9D = {z G C : = 1}. One adds a pure point £ £ dD to v to form the new probability measure 

u = (i_ 7 ) z/ + 7 ^ j l> 7 >0. (3.3) 

Then one has the following classic result relating the original and the perturbed orthogonal poly- 
nomials, which is an analog of (|2.4p in Theorem [ 



Theorem 3.1. (Geronimus [HJdSJJ Suppose the probability measure v is defined as in 113. 3\) . Then 
the n-th monic orthogonal polynomial of v is given by 

where & n (z) and (p n (z) are the monic and normalized orthogonal polynomials of u respectively; and 



z 



The point mass problem on the unit circle was further investigated by Cachafeiro-Marcellan [2] 
and Simon (see Chapter 10.13 of [21]) using very different approaches. In |27| Wong applied the 
Christoffel-Darboux formula to Simon's result and proved the point mass formula, which shows that 
the recurrence coefficients (a n )'^' =0 of du and dv are related as a n (du) = a n (dv) + A n ((), where 



A fA\ (1 ~ \an\ 2 ) 1/2 (p n +l(0<Pn(Q v 1 1 fs ST I ( t\V /q r\ 



The reader may compare (|3,5p with the formulae for a n and b n in Theorem! 

In [SHJEHIES] Wong studied the point mass problem on the unit circle. In particular, the class W p 
consisting of measures of which the Verblunsky coefficients are of p-generalized bounded variation 
was identified. It was proven that upon adding a point mass to a measure in W p , we obtain a new 
probability measure in W p +i. Inductively, we can add a finite number k of distinct pure points to 
a measure in W p one after another and we will end up with a new probability measure in W p +k- 
These results on the unit circle could be read in parallel to Theorem [ 
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4. The Transfer Matrix Approach 

4.1. Structure of this section. In Section ^. 2l we consider the recurrence relation in matrix form. 
Since xo is not in the support of the measure, Aj(xo) is hyperbolic for all large j. Because of the 
hyperbolicity and the fact that Aj — > Aao, we can apply Kooman's Theorem to prove that there is 
an analytic function [Zona neighborhood of A^ such that for all large j, Aj can be simultaneously 
diagonalized (see Section [43]) . 

In Section 14. 4| we consider a representation that stems out from the diagonalization. In Propo- 
sition I4.2p , we show a trichotomy based on that representation and prove Theorem 12.11 for two of 
the three cases, which are special cases that are easier to handle. 

The remaining case, being the most difficult one, will be treated in Section T4.51 We prove several 
estimates for this particular case. In Section [4. 6\ we combine those estimates to obtain asymptotic 
formulae for p n (xo). 

A diagram summarizing the results is provided at the end of the section for the convenience of 
the reader. 



From now on, let fi be a measure in Ai(b, a) with recurrence coefficients satisfying 
(I2TT0D and fl2TT|) . 

4.2. The recurrence relation and p n (xo). Note that under the conditions that a n — > a ^ and 
b n -> b, 

lim A j (x) = A 00 (x):=a- 1 ( X - 2 b . (4.1) 



A ao (x) has eigenvalues 



± (x-b)± y/(x - b) 2 - 4a 2 



2- ' < 42 » 

Therefore, A^x) has distinct eigenvalues in R if and only if x £ R\[6 — 2a, b + 2a]. In that case, 
one of which has absolute value strictly greater than 1 and the other strictly less than 1. We say 
that A^x) is hyperbolic. 

Now consider a fixed point xq G R\[6 — 2a, b + 2a]. Without loss of generality, we assume 
x > b + 2a. By P~2|) , we have 

A+ > 1 > A" (4.3) 

for all large n. 

From now on we will write Aj(xo) as Aj and similarly for other objects that appear in the proof. 

4.3. Hyperbolicity of Aj and Kooman's Theorem. Since Aj —> Aoo, Aj is also hyperbolic for 
all large j. This allows us to use the following result by Kooman. Adapted to suit the context of 
this paper, the theorem reads as follows: 

Theorem 4.1 (Kooman |13} I14j). Let A be an t x t matrix with distinct eigenvalues. Then there 
exist e > and analytic functions U{B) and D(B) defined on S e = {B : \\B — A\\ < e} such that 

(1) B = UbDbU^ 1 , Db commutes with A. 

(2) Ub is invertible for all B G S e . 

(3) U A = 1,D A = A. 

(4) By picking a basis such that A is diagonal, we can have all Db diagonal with entries being the 
eigenvalues of B. 

Remarks: 

(1) The formulation of Theorem \4-l\ is similar to Theorem 12.1.7 of [21] . except that in |21j the 
statement was meant for quasi-unitary matrices. In fact, the proof holds as long as A has 
distinct eigenvalues. 
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(2) Kooman's Theorem first appeared in Theorem 1.3 o/[13]. The first application of Kooman's 
Theorem to orthogonal polynomials was made by Golinskii-Nevai |10j for the unit circle 
case. They proved that if the recurrence coefficients of the measure fi on <9B satisfy a n — > 
and ifY^n ll^-n+l — < 00 > ^ en th> e a - c - P ar t of the measure is positive almost everywhere 
on <9B. 



Let G be the matrix that diagonalizes Aoo, i.e., 

A x = G^DG, (4.4) 

where 

A+ 



D ={o X-) ^ 

and X^ are the eigenvalues of A.^ (defined in (|4.2p ). Pick a basis in which j4oo is diagonal. Then 
by the construction of the function D in Theorem 14.11 above, there exists an integer 

N > N(e) (4.6) 

such that Aj is in some S e neighborhood of A^ and Da. is a diagonal hyperbolic matrix under this 
basis for all j > N. In other words, there exist diagonal matrices 

such that Da^xo) = GDjG~ x and the eigenvalues A^ satisfy 



In fact, by (|1.12p and a straightforward computation, we can show that A^- are roots of the charao 



A+ > 1 > AT, lim Xf = A ± . (4.8) 

id a straightfor 

teristic polynomial 

z 2 - (s -bj)z + a) = 0. (4.9) 

Next, we define 

Gj = U Aj G. (4.10) 
Then for j > N, the matrix Aj can be expressed as 

A j = G j D j Gj\ (4.11) 

4.4. Representation of the recurrence relation. Let E be an integer, we will choose E later 
in the proof (see Proposition I4.2p . 

For n > E > N , the transfer matrix applied to the basis vector gives 



T n — G n D n G n 1 G n -iD n -iG n 1 : 1 ■ ■ ■ De+\G e \. x GeVe, 



(4.12) 



where 



ve = ( V \ E) ) := DeG^Ae-i ••• A x ( ) . ( 1. 13) 



Following (I4.12p . we consider the equation 



D n G- n x G n ^D n ^G-\ ■ ■ ■ De+iG^Gev = L n UnV \ E) ] , (4.14) 

\WnV y 2 ' 
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where u n and w n are defined implicitly by ()4.14p above and 

n 

L n= II A+ for n>E. (4.15) 

k=E+l 

Also, we define 

( V t)) = v(n) = L.H.S. of gUD = L n ( UnV \^ ] . (4.16) 

Proposition 4.1. For any E > N, either v\ or v\ in (|4.13p is non-zero. 
Proof. Observe that by (14.131) . 



Since Gjv is invertible and that a n > for all n, t; = implies pn{xq) = pn-i(xo) = 0. This 
contradicts the fact that the zeros of p n (x) and p n _i(x) strictly interlace for all n and 

□ 

In the following proposition, we are going to identify the trichotomy about the pair m 
(I4THD : 

Proposition 4.2. Lei uj; and fr, ^ e defined as in (14.130 and iV 6e defined in (I4.6|) . T/ien one 
o/ i/ie following is true: 

( E) ( E) 

(1) For some E > N , both v\ and v\ are non-zero. 

(2) For all E> N, v { 2 E) = 0. In that case, 



\k=N / 



Pn(x ). (4.18) 



(3) For all E > N, v[ E) = 0. In that case, 

/i-i 



Pj(xo) = J] X k) M*o). (4.19) 



Kk=N 



Remark: Case (1) will be treated in Section \4-5\ below. 

Proof. Suppose for all E > N , v\ =0. By Proposition 14. 1\ v\ ^ for all E. That is equivalent 



to 



^ f Pe(x ) 



Ge v i ) = ( PE ^\ A, VE>N. (4.20) 
V J \a E -iPE-i(x ) J y ' 

Recall from (14.111) that Ae = GeDeG^ 1 ■ ^ n °f ner words, Ge is the change of basis matrix that 
maps the vector (1,0) to the eigenvector of Ae with eigenvalue A^,. Therefore, 

PE + i(xo)\ =A J Pe(x ) ) =x+ ( Pe(xo) V (421) 
aEPE{x )J \a E -iPE-i{x ) J * \a E -iPE~i(xo) J 

By an inductive argument, we obtain (|4.18p . 

The proof for (3) is identical in nature except that we have 

(pe + i(x )\ f pe(x ) \ (422) 

\a E PE{xo)J a \a E -iPE-i{xo) J 

instead of (14.21H because Ge maps (0, 1) to the eigenvector of Ae with eigenvalue A^. This proves 
flUS}. 
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□ 

Now that we have asymptotic formulae for p n (xo) for Cases (2) and (3), we are going to consider 
Case (1). In fact, if there exists E such that v\ and are non-zero, the asymptotic formula for 
Pn(xo) will be quite similar to (|4.18p or (I4.19h except for the presence of error terms, which will be 
analyzed in the next section. 



4.5. Several estimates. By Proposition 14.21 there are three cases. In this section we focus on 
Case (3): i.e., there exists an E > N such that v[ and v\ are both non-zero. 



Consider a fixed E. For convenience, we shall write 



(E) (E) , 

v\ = v\ ,V2 = v 2 and v = ve- 



Proposition 4.3. There is a constant C such that 

\\v{n + 1) - D n+1 v(n)\\ < C\\A n+ i - A n \\\L n \ (\u n \ + \w n \) . 

Proof. Note that 

v(n + 1) - D n+1 v(n) = D n+1 (G~ +1 G n - l) v{n). 

The goal is to give bounds for each of the components on the right hand side of ()4.25p . Since U 
is analytic on S e , on some compact subset of S e there exist constants 771, 772 > such that 



(4.23) 

(4.24) 
(4.25) 



and 

Thus, for 7] = 771772, 
Furthermore, 
and 



||G n+ i - G n \\ < ||G||||C/An+i ~ U A„\\ < m\\An+l - An\\ 

l|G;i 1 ||<i|G- 1 ni|c/- t ; +i ||< 77 2 . 

G n ^i 1 G n — 1|| = ||G n +i (G n — G n +i) || < 77||A n+ i — A n \ 



sup \\D n \\ = sup |A„ I < 2|A 

n>N n>N 



vinj 



u n L n v 1 
w n L n v 2 



< Ci\L n \ (\u n \ + \w n \ 



(4.26) 
(4.27) 

(4.28) 
(4.29) 

(4.30) 



where Ci = max{ \ vi |, \v 2 \} > 0. By applying (EQgj) . (11391) and (jPUj) to (14T251) . we finish the proof 
of Proposition 14.31 □ 



Proposition 4.4. Let u n and w n be defined as in (|4.14p and v\,v% i= 0. Then the following 
inequalities hold: 

(1) There is a constant C3 such that 

\u n+1 - u n \ < C 3 ||^ n+ i - A n \\ (\u n \ + . (4.31) 

(2) There is a constant C4 such that 

a; 



w n+ i 



A 



' n+l 
+ 

n+l 



W, 



< C 4 \\A n+ i - A n \\ (\u n \ + \w n \) . 



(4.32) 
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Proof. Recall that L n+ \ = X^ +1 L n and Vn = L n u n v\. Since v\ 7^ 0, 



\Un+l ~ Ur. 



n 1 - A+ v v 

J n+1 A n+l u n 



(1) 



V\L. 



n+l 



< 



\w(n + 1) — D n+ iw(n) 
\v x L n+ i\ 



(4.33) 



(2) 

Similarly, since v n = L n w n V2 and v<i ^ 0, 
A 



n+l 
- — '"n 



A 



n+l 



J 2 ) _A" v {2) 

u n+l A n+l u n 



< 



\w(n + 1) - D n+ iw(n) 
\v2L n +i\ 



Apply Proposition 14.31 to the equations above to obtain (|4.3ip and (|4.32[) . 



(4.34) 



□ 



The following lemma concerning u n and w n is central for this paper. As we shall see in the proof 
of Theorem 12.31 it implies the dichotomy between exponential decay and exponential growth of 

Pn(xo)- 

Lemma 4.1. Let u n and w n be defined as in (14. 14ft and vi,i)2 ^ 0. Then one of the following is 
true: 

(1) There exists a constant C such that \u n \ < C\w n \. Moreover, given any e > 0, there exist 
an integer N e and a constant C e such that 



\Wn\ < C e 

(2) \w n /u n \ — > 0. Furthermore, Uqq = 



A+ 



+ e 



Vn > N r . 



(4.35) 



liim 



u n exists and it is non-zero. 



Proof. There are two possible situations concerning u n and w n : 

(1) There exist a fixed integer K$ and a constant C such that |it n | < C|io n | for all n > Kq. 

(2) For any integer K and any constant H, there exists an integer tik,h > K such that \u nK B \ > 



H\w. 



Case (1): By ()4.32p . for n > max{N, Kq}, there is a positive constant Cj such that 

A, 



\Wn+l 



< 



An 



+ CV||j4 n +i — A r 



W r , 



(4.36) 



Recall that ||^4 ra +i — A n \\ —> and — > A 1 * 1 . Thus, given any e > 0, there exist an integer N e 
and a constant C f such that 



\w n \ < C f 



A+ 



+ e 



Vn > N e . 



(4.37) 



In other words, w n decays exponentially fast. Hence, u n also decays exponentially fast to zero given 
the inequality \u n \ < C|w n |. This corresponds to (2a) of Lemma 14. 11 



Case (2): Let r n = w n /u n . The statement r n — > is by definition equivalent to proving that given 
any e > there exists an integer J e such that \rj\ < e for all j > J t . 

First, we show that both u n and u n+ \ are non-zero, as (|4.39p below will involve u n and u n+ \ in 
the denominator. 

By the assumption we can choose any H, we choose one such that 1/H < e. Consider any 
fixed pair {K,H) (the choice of K will be made later in the proof). The existence of an integer 
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n = n,K,H > K such that \r n \ < 1/H = e is guaranteed, which implies u n ^ 0. Furthermore, by the 
triangle inequality and (|4.3ip . 



Un+l 



U,< 



> 1 



u n +i - u n 



Un 



(4.38) 



> l-C 3 \\A n+1 -A„||(l + |r n |) > 



which implies that u n +i is also non-zero. 
Next, observe that 



< 



r n +i ~ — fr Tl 

An 

Wn+l K w n 



u n+l A n U n+ i 

Wn+i - (A"/A+)u; 



+ 



Wn 



U n +1 



+ 



K_ 



u n +i u n 

U n ~ U n+ i 



(4.39) 



U n +l 



By (I4.3ip and (I4.32p . there is a positive constant Cs such that 

A n 

An 

< 1 + l^l'A- /A+ 1 c4An+i _ An \\ {lUn \ + K[) 



(4.40) 



\u n +i\ 

C 8 (l + |r„||A-/A+|)p„+i- A r , 



I Wn+l | 



By inverting (|4.38p . 



< 



1 



Un+l 

Then we plug it into ([4.401) t° 

A„ 



Fn+l| < 



+ 



1 - C 3 \\A n+l - A n \\(l + \r n 
C 8 (l + |r n ||A-/A+|)(l + |r ri |) 



A 



n+l 



A, 



(4.41) 



(4.42) 



l-C 3 \\A n+l -A n ||(1 + \r n \ 

Since ||A n +i — A n \\ — > 0, the second term on the right hand side of (j4.42p can be arbitrarily 
small if n is sufficiently large. Hence, for any sufficiently large K, there exists n > K such that 
|r n+ i| < \r n \ < e. 

Inductively, we can aply the same argument to r n+ \ to prove that |r n+ 2| < e - As a result, \rj\ < e 
for all large j. This proves \w n /u n \ — > 0. 



Next, we are going to show that limn^oo u n exists. Divide (|4.3ip by \u n \. Since \r n \ — > 0, 

u n +i 



1 



U n 



< C\\A n+l - A n \\ (1 + \r r , 



0. 



(4.43) 



Moreover, since log is analytic near 1, in an e-neighborhood of 1 there is a constant E such that 

| \ogz\ = | log C - log 1| < E\z - 1|. (4.44) 

By (0331), 



log 



U n +1 



U, 



<c||Ah-i-A»I|. 



(4.45) 
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Therefore, the series ^^jL^log (uj + i/uj) is absolutely convergent. Furthermore, as we have seen 
in ()4.38p . Uj 7^ for all large j. Thus, log Uj is well-defined and the following limit 



71-1 



lim logu n = lim } (logtu+i — logUj) +log 



3=P 



exists and is finite. 

This corresponds to the second part of (2b) and concludes the proof of Lemma 14.11 



(4.46) 



□ 



4.6. Implications of Lemma 14.11 By Lemma 14.11 there are two possible situations: 

Case (1). Suppose (2a) of Lemma l4.1l is true. We will prove that this corresponds to the case 
/j,(x ) > 0. 

First, observe that 



Pn{X0 \ xW B (z )ft 
a n -iPn-i(x ) J ny u; 10 



G n L n 



U n Vl 

w n v 2 



Hence, given any e > 0, there exists a constant K t > such that 



(4.47) 



T n (x ) 



< \\G r 



n w 

k=N+l 



U n Vl 

w n v 2 



<KA\\ + \+e) 



+ e 



Since |A | < 1 and |A /A + | < 1, p n (xo) goes to zero exponentially fast by (|4,48|) . 
It is well-known (see for example [5]) that 



N 



-1 



lim Vp n (3;o) 

N^roo ' — ' 



ra=0 



(4.48) 



(4.49) 



Hence, n(xo) > and we are just varying the weight of xq. It is easy to prove that by (|2.8p and 
(|2.9p . both a n — a n and b n+ \ — b n+ \ go to exponentially fast. 

This is in agreement with Simon's result (Corollary 24.4 of [22J) that varying the weight of a pure 
point will result in exponentially small perturbation of the recurrence coefficients. 



Case (2). Suppose (2b) of Lemma 14. II is true. Let 

Gn= (9l,n 9' } , n \ G= (9l £\ (45Q) 
\92,n 92,J \92 9 2 J ^ ' 

Since p n (xo) is the first component of the vector G n L n (u n vi, w n V2) T , we have 

p n (x ) = L n (g 1:n u n vi + g' ln w n v 2 ) 

= L n u n (gi tn + g' ltn r n v 2 ) (4.51) 

= Ln{Uoo9lVl + O(l)). 

The last equality holds because r n = u n /w n — > by (2b) of Lemma liTTl 
Here is a summary of the results in this section: 
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vi = 





Definition 

of (vi,V 2) 








Vi,V 2 / 







Lemma 14.11 




5. Proof of Theorem 12.11 

By the discussion in Section^ p n {%o) 2 is either exponentially increasing or exponentially decaying 
towards zero. Moreover, by (|4.49|) above, 

oo 

fl{x ) > ^2pn{xo) 2 < OO. (5.1) 

n=l 

Therefore, fJ,(xo) > if and only if p n (xo) 2 is exponentially decaying towards zero, which cor- 
responds to ()4.19j) and (|4.48|) ; fJ,(xo) = if and only if p n (xo) 2 is exponentially increasing, which 
corresponds to (|4.18|) and (|4.5ip . 



6. Proof of Theorem 12.21 
Proof. Let K n (x,y) be the reproducing kernel of the measure p, which is given by 

n 

K n (x,y) = ^2p j (x)p j (y). (6.1) 
3=0 

Since P n (x) is a polynomial of degree n, 

P n (x) = I P n (y)K n (x,y)dp(y) 

(6-2) 

P n (y)K n (x, y)dp(y) - -yP n (x )K n (x, x ) . 
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Moreover, P n (y) is orthogonal to all polynomials with degree < n — 1 with respect to the inner 
product ( , )du- Therefore, 

J Pn(y)K n (x,y)dfl(y) = j P n (y)p n (y)p n (x)dfi(y) 

= Pn(x) J P n {y)p n (y)dil{y) 

and that 



(6.3) 



(P n (y),p n (y)) = (P n (y), K, n y n ) = . (6.4) 

Therefore, by (|6.2p and (|6.3D . 



- 2 



Pn{x) = ^ Pn(x) ~ jP n (x )K n (x,X ). (6.5) 

Now plug in x = xq into (16. 5^ . Upon rearranging, we get 

\Hi n J l + jK n (x ,x ) 

Putting ()6.6p into (|6.5p . we arrive at f|2.4|) . In particular, note that both P n (x) and P n (x) are 
monic polynomials. Hence, by comparing the coefficients of x n on each side of (12.41) . we get 

i = f^Vfi- 1Pn{x ; )2 (6.7) 

\K n / V 1 + 7A„(xo,x )/ 

«nV 1 + 7-^n-l (^0,^0 ) in q\ 

— -;- — s-- (6.8) 

K n J l + jK n (xo,x ) 

Recall that a n = n n -i/K n . This proves (|2.7j) . 

Next, we prove the formula for b n +\. Let m n be the coefficient of x n ~ l in P n (x), i.e., 

P n (a;) = x n + m n x n ~ l + lower order terms . (6-9) 

By the recurrence relation (11. 6p . 6 n +i is given by the coefficient of x n in xP n (x) — P n +i(x), which 
can also be expressed as 

b n+1 = m n - m n+1 . (6.10) 
To prove formula (I2.9P for 6 n +i, we will compute m n . By (j2.4|> . 

rtin = ( ^ ] (m n - 7jn_3fr_ [p n ( Xo ) Knmjl +p n _i(j; )K ri _i] J . (6.11) 

VK n y V 1 + jK n (x ,x ) ) 

The coefficients of m n in (|6.1ip are given by 

— V ( i - ^"( x q) 2 A / 6 12 ) 



which is equal to 1 by (j6.7j) . Therefore, 



7P n (x )p n -l(^o) K n— 1 



1 + 7^T n (xo,xo) 

7P n (x )p„-i(x ) 
1 + 7i^ n _i(j;o,so) 



(6.13) 



The last equality follows from the expression of (ft n /K n ) 2 in (16. 8p . 

This concludes the proof of Theorem 12.21 □ 
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7. Proof of Theorem 12.21 
We are going to separate the proof into two different cases. 

Case (1): xq is a pure point of \x. This is the easier case. Since 

lim K n (x ,x ) = /i(_o) -1 , (7.1) 

n— >oo 

it is clear by (|2.8p that t n — > 1 and a n — > a as n — > oo. Furthermore, recall that P n+ \(xo) = 
p n+ i(x ) / n n+1 and a n+1 = n n /n n+1 . Hence, 

h , = jPn+i(x )p n {x )K n _ ja n+1 p n+1 (xo)p n (x ) ^ Q 
l + jK n (x ,x ) 1 + -fK n (x ,x ) 

because p n (xo) —> exponentially fast. Therefore, 

b n+ i = b n+ i - h n _i + h n -)• 6 as n -)■ oo. (7.3) 
It is trivial to prove (I2.1ip for this particular case so we shall omit the proof. 



Case (2): xq is not a pure point of /i. To prove that a n — > a, we are going to prove that 
lim^^oo t n exists by employing the following theorem: 

Theorem 7.1 (Cesaro-Stolz [3]). Let (r n ) nG ^, (K n ) n£ f$ be two sequences of numbers such that K n 
is strictly increasing and tends to infinity. If the following limit 

lim Tn ~ ^ (7.4) 

exists then it is equal to linin._j.oo ^n/Kn- 

Let T n = 1 + ~/K n _i(x , x ) and K n = 1 + ^/K n (xQ, xq). Since A+ > 1, L n — > oo. Hence, 
Pn(^o) - > oo, which also implies that lf n is an increasing sequence to infinity. This allows us to 
apply the Cesaro-Stolz Theorem to prove that 

,. , ,. r n -r n _i p n _i(x ) 2 i 

hm t n = hm — = lim — — - = — 5 - (7.5) 

The last equality of (|7.5p follows from the asymptotic formula (|4.5ip for p n (xo)- As a result, 
lim ra _ >00 a n = lim r) ,_s. 00 a n = a. 



Next, we are going to prove that lim^oo b n = b. We shall employ the Cesaro-Stolz Theorem 
again. Let 

T n = IP n+i(xo)p n (x ) and K n = 1 + jK n (x ,x ). (7.6) 
By the asymptotic formula (I4.5ip for p n (xo), 

F n - T n _i = 7Pn(^o) (gn+l(go) -Pn-l(x )) + _ J_ , > 

^n-^n-1 7_n(x ) 2 A+ ' l ' j 

Therefore, 

lim h n = lim ^ + iPn + i(x )p n (x ) = J X+ _1\ m 

n.-5-oo n-5-oo I +jK n (x ,Xo) \ A+ / 

and by (|7.3p . 6 n — > b. 

We are going to prove that the sequence (a n ) n is of bounded variation. 
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Since liuin^^ a n = a / 0, for all large n there exists a constant C such that 
I 2 



a 2 n+ i — a n | — \a n+ \ — a n \\a n+ i + a n \ > C\a n+ i — a n \. (7.9) 



Therefore, it is enough to prove that (a n ) n is of bounded variation. Furthermore, by the formula 
for a n (|2.7p and the fact that lim^oo t n = 1/A + , it suffices to show that the sequence (t n ) n is of 
bounded variation. 

To do that, we are going to show that the sequence (r n ) n is of bounded variation, where 

Pn(xp) 2 , s 

T n = _! , 77 , r • (7.10) 

7 + K n {xo,x ) 

There are two cases to consider, p n (xo) being in the form (|4.18|) or p n (xo) being in the form 
(jSU). 

If p n (xo) is in the form ()4.18j) . 

Ll Pn{x ) 2 



7- 1 + K n (x ,x ) \+ 2 

If p n (xo) is in the form (|4.5ip . 



(7.11) 



L 



n 



n TnV2) ■ (7.12) 



7 1 + K n (xo,x ) s_ 



5 

Observe that Pn{xq) 2 /(X^) 2 is of bounded variation, because pn{xq) 2 isaconstant and (1/(A+) 



2^ 



is of bounded variation. Moreover, |<7i,n+i — 9i,n\, \g'i n +i ~ 9*2 n\ an d l r n+i — f n \ are of the order 
0(||G n+ i — G n ||), which implies they are also of bounded variation. Therefore, given (|7,lip and 
()7.12p . it remains to show that L 2 l /('y~ 1 + K n (xo, xo)) is of bounded variation. 
We will make use of the simple fact: given the equality 

— - — = y^-y^ (7.13) 

Vn Vn+1 Vn+lVn 

if lim n _ i . 00 y n = y 7^ 0, (y n ) n is of bounded variation if and only if (l/y n ) n is of bounded variation. 

Hence, we will prove that (7 _1 + K n (xo, xq)) / X n is of bounded variation and its limit exist when 
n goes to infinity. To prove the latter, observe that by the Cesaro-Stolz Theorem, 

r 2 - r 2 t 2 ((\ + \ 2 - \\ 

lira K+ \ L » r = K UA " + , lJ lj = t ((A+) 2 - 1) * 0, (7.14) 

n^oo K n+1 (x ,Xo) ~ K n (x ,X ) p n+ i(X y 

where £ is a non-zero constant whose value depends on whether p n (xo) is in the form (|4.18|) or 
(I43TT) . 

For the convenience of computation we will define a few more objects below. First, we let 

if " £Ar + 1 . (7.15) 
if < n < N v ' 

Then by (I4.15p . L n = 11^=1 Ar Moreover, recall the definition of u n in (I4.14p . which was only 

defined for n > N. For < n < N, let u n w n be defined implicitly by (I4.14|) . The introduction of 
these objects will not affect the result of our computation. 
Observe that 

7' 1 + K n (x ,x ) 1 + 7' 1 , p / 71R n 
ft = — 72 — + R n, (7.16) 

where 

R n =7TEftW = E^^+ £r . (7-17) 

j=1 j=1 A j+1 • • • An 
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with the convention that A^ +1 • • • A+ = 1 when j > n. 
Let 



K n _l{XQ,XQ) 

L n-1 



^ ^(ffi^i+ff^r^) 2 

j=0 A j+1 A n-1 



(7.18) 



Then 

7 _1 + K n (x ,x ) 7 _1 + K n _i(x ,xo) 



LI 



< a 2( V +7 'L +l^n-gnl- (7.19) 

(A n _i^ • • • A ^r 



Recall that A+ — > A + > 1. Hence, 



E 



1 



(A„ + ---A + ) i 



few 



< oo. 



(7.20) 



Thus, the first sum on the right hand side of (|7.19p is summable. 

Next, observe that upon rearranging the indices of S n in (|7.18p . we have 



Rn >Sn | 

^2 u ^ 9l 'i Vl + 9i,j r i V2 ) 2 u2 j-i(gi,j-i v i + a'i d -ifj-iV2) 2 



(a;---a+_i) 2 



n 2| _ 
\ ^ u j\ e j ~ e J — 1 1 



E 



1 l„,2 ,,2 



(A+---A+)2 (A+---A+_i) 5 



(7.21) 



where 



ej = (52j^i + g' 2j rjV 2 f 



(7.22) 



Now we proceed to estimate each of the sums on the last line of (|7.2ip . Recall that \uj — Uj—i\ = 
0(\\Gj — Gj-i\\). Moreover, Uj — > u^. Hence, YE- Wj ~ u j-i\ < 00 • Therefore, for some constant 
C, 



X " "] "j , ', i 



XT / \ + \ + 

n=lf=l ^ i A «-l' ) 



sc EK 2 -«l-il Etitw <«-■ 



Similarly, 



because 



00 n „,2 



X/X/ c \ . . . \+ r: 

n=l 7 =1 V A i+l A ™ i 



< OO 



ki - ej-il = 0(|r, - r,-_i|) + 0(\g hj - + 0(1^ - tfij-xl) 

= 0(||G i -G,_ 1 ||). 

Finally, we consider the second sum on the right hand side of (I7.2ip . Note that 

n-l 

|Af - Af | < C\\f - A+| < |A+ +1 - A+| , 



(7.23) 



(7.24) 



(7.25) 



(7.26) 
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where C is a positive constant independent of j and n. As a result, 

1 1 



EE 

n=l j=l 

oo n 



< 



1 



(7.27) 



oo n 



n-1 



C EE 



X] l A fe+i ~ A * 



n=l j=l ( A i+1 • • • An fc=j 

Consider a fixed fcsN. We count the coefficients of |A^ +1 — in the last line of (|7,27p . From 
the equation, we know that j < k < n. Thus, the coefficient is 

oo k 

^ (\+ \ + \2 

n =k+l j=l 



3=1 \n=k+l ( A j+1 ---A ») 2 < 



(7.28) 



E 



E 



which is bounded above by a constant B independent of k. Therefore, 

oo 

dZZZD < Ci? ^ | A+ +1 - A+ 1 < oo. 



(7.29) 



k=l 



Going back to (j7.21|) . we conclude that Y^ n \ R n — S n \ < oo. Hence, a n is of bounded variation. 

Finally, we are going to show that b n is of bounded variation. Since a n is of bounded variation, 
by (|2.9p . it suffices to show that (h n ) n is of bounded variation, h n as defined in (|7.2f) . 

By 631]), 



h r , 



[gi,n+ivi + g' 1 , n+1 r n+ iv 2 ) (gi, n vi +g'i :n r n v 2 ) ■ 



(7.30) 



7 1 + K n (x ,x ) 

Thus, it boils down to proving that L 2 n j (7 _1 + K u (xq, xq)) is of bounded variation. To that end, 
we use the same argument following (|7.10p . 
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